A Novel Mutation in the INSR Gene Causes Severe Insulin Resistance and Rabson-Mendenhall Syndrome in a Paraguayan Patient.

Rabson-Mendenhall syndrome (RMS) is a rare autosomal recessive disorder characterized by severe insulin resistance, resulting in early-onset diabetes mellitus. We report the first case of RMS in a Paraguayan patient. The patient is a 6-year-old girl who presented with hypertrichosis, acanthosis nigricans, nephrocalcinosis, and elevated levels of glucose and insulin that served as diagnostic indicators for RMS. Genetic testing by next-generation sequencing (NGS) revealed two pathogenic variants in exons 2 and 19 of the INSR gene: c.332G>T (p.Gly111Val) and c.3485C>T (p.Ala1162Val), in combined heterozygosis. The novel INSR c. 332G>T variant leads to the substitution of glycine to valine at position 111 in the protein, and multiple in silico software programs predicted it as pathogenic. The c.3485C>T variant leads to the substitution of alanine to valine at position 1162 in the protein previously described for insulin resistance and RMS. The management of RMS is particularly challenging in children, and the use of metformin is often limited by its side effects. The patient was managed with nutritional measures due to the early age of onset. This report expands the knowledge of RMS to the Paraguayan population and adds a novel pathogenic variant to the existing literature.

Effect of Essential Oil Components on the Activity of Steroidogenic Cytochrome P450.

Endocrine-disrupting chemicals (EDCs) may impact the development of prostate cancer (PCa) by altering the steroid metabolism. Although their exact mechanism of action in controlling tumor growth is not known, EDCs may inhibit steroidogenic enzymes such as CYP17A1 or CYP19A1 which are involved in the production of androgens or estrogens. High levels of circulating androgens are linked to PCa in men and Polycystic Ovary Syndrome (PCOS) in women. Essential oils or their metabolites, like lavender oil and tea tree oil, have been reported to act as potential EDCs and contribute towards sex steroid imbalance in cases of prepubertal gynecomastia in boys and premature thelarche in girls due to the exposure to lavender-based fragrances. We screened a range of EO components to determine their effects on CYP17A1 and CYP19A1. Computational docking was performed to predict the binding of essential oils with CYP17A1 and CYP19A1. Functional assays were performed using the radiolabeled substrates or Liquid Chromatography-High-Resolution Mass Spectrometry and cell viability assays were carried out in LNCaP cells. Many of the tested compounds bind close to the active site of CYP17A1, and (+)-Cedrol had the best binding with CYP17A1 and CYP19A1. Eucalyptol, Dihydro-ß-Ionone, and (-)-α-pinene showed 20% to 40% inhibition of dehydroepiandrosterone production; and some compounds also effected CYP19A1. Extensive use of these essential oils in various beauty and hygiene products is common, but only limited knowledge about their potential detrimental side effects exists. Our results suggest that prolonged exposure to some of these essential oils may result in steroid imbalances. On the other hand, due to their effect on lowering androgen output and ability to bind at the active site of steroidogenic cytochrome P450s, these compounds may provide design ideas for novel compounds against hyperandrogenic disorders such as PCa and PCOS.

Pharmacogenomics of poor drug metabolism in greyhounds: Canine P450 oxidoreductase genetic variation, breed heterogeneity, and functional characterization.

Greyhounds metabolize cytochrome P450 (CYP) 2B11 substrates more slowly than other dog breeds. However, CYP2B11 gene variants associated with decreased CYP2B11 expression do not fully explain reduced CYP2B11 activity in this breed. P450 oxidoreductase (POR) is an essential redox partner for all CYPs. POR protein variants can enhance or repress CYP enzyme function in a CYP isoform and substrate dependent manner. The study objectives were to identify POR protein variants in greyhounds and determine their effect on coexpressed CYP2B11 and CYP2D15 enzyme function. Gene sequencing identified two missense variants (Glu315Gln and Asp570Glu) forming four alleles, POR-H1 (reference), POR-H2 (570Glu), POR-H3 (315Gln, 570Glu) and POR-H4 (315Gln). Out of 68 dog breeds surveyed, POR-H2 was widely distributed across multiple breeds, while POR-H3 was largely restricted to greyhounds and Scottish deerhounds (35% allele frequencies), and POR-H4 was rare. Three-dimensional protein structure modelling indicated significant effects of Glu315Gln (but not Asp570Glu) on protein flexibility through loss of a salt bridge between Glu315 and Arg519. Recombinant POR-H1 (reference) and each POR variant (H2-H4) were expressed alone or with CYP2B11 or CYP2D15 in insect cells. No substantial effects on POR protein expression or enzyme activity (cytochrome c reduction) were observed for any POR variant (versus POR-H1) when expressed alone or with CYP2B11 or CYP2D15. Furthermore, there were no effects on CYP2B11 or CYP2D15 protein expression, or on CYP2D15 enzyme kinetics by any POR variant (versus POR-H1). However, Vmax values for 7-benzyloxyresorufin, propofol and bupropion oxidation by CYP2B11 were significantly reduced by coexpression with POR-H3 (by 34-37%) and POR-H4 (by 65-72%) compared with POR-H1. Km values were unaffected. Our results indicate that the Glu315Gln mutation (common to POR-H3 and POR-H4) reduces CYP2B11 enzyme function without affecting at least one other major canine hepatic P450 (CYP2D15). Additional in vivo studies are warranted to confirm these findings.

Exploring Novel Variants of the Cytochrome P450 Reductase Gene (POR) from the Genome Aggregation Database by Integrating Bioinformatic Tools and Functional Assays.

Cytochrome P450 oxidoreductase (POR) is an essential redox partner for steroid and drug-metabolizing cytochromes P450 located in the endoplasmic reticulum. Mutations in POR lead to metabolic disorders, including congenital adrenal hyperplasia, and affect the metabolism of steroids, drugs, and xenobiotics. In this study, we examined approximately 450 missense variants of the POR gene listed in the Genome Aggregation Database (gnomAD) using eleven different in silico prediction tools. We found that 64 novel variants were consistently predicted to be disease-causing by most tools. To validate our findings, we conducted a population analysis and selected two variations in POR for further investigation. The human POR wild type and the R268W and L577P variants were expressed in bacteria and subjected to enzyme kinetic assays using a model substrate. We also examined the activities of several cytochrome P450 proteins in the presence of POR (WT or variants) by combining P450 and reductase proteins in liposomes. We observed a decrease in enzymatic activities (ranging from 35% to 85%) of key drug-metabolizing enzymes, supported by POR variants R288W and L577P compared to WT-POR. These results validate our approach of curating a vast amount of data from genome projects and provide an updated and reliable reference for diagnosing POR deficiency.

Exploring the Potential of Sulfur Moieties in Compounds Inhibiting Steroidogenesis.

This study reports on the synthesis and evaluation of novel compounds replacing the nitrogen-containing heterocyclic ring on the chemical backbone structure of cytochrome P450 17α-hydroxylase/12,20-lyase (CYP17A1) inhibitors with a phenyl bearing a sulfur-based substituent. Initial screening revealed compounds with marked inhibition of CYP17A1 activity. The selectivity of compounds was thereafter determined against cytochrome P450 21-hydroxylase, cytochrome P450 3A4, and cytochrome P450 oxidoreductase. Additionally, the compounds showed weak inhibitory activity against aldo-keto reductase 1C3 (AKR1C3). The compounds' impact on steroid hormone levels was also assessed, with some notable modulatory effects observed. This work paves the way for developing more potent dual inhibitors specifically targeting CYP17A1 and AKR1C3.

Expanding the p.(Arg85Trp) Variant-Specific Phenotype of HNF4A: Features of Glycogen Storage Disease, Liver Cirrhosis, Impaired Mitochondrial Function, and Glomerular Changes.

Introduction: The p.(Arg85Trp) variant-specific phenotype of hepatocyte nuclear factor 4 alpha shows a complex clinical picture affecting three different organ systems and their corresponding metabolisms. Little is known about the molecular mechanisms involved and their relationship with the diverse symptoms seen in the context of this specific variant. Here, we present data of a new patient that expand the clinical phenotype, suggesting possible disease mechanisms. Case Presentation: Clinical data were extracted from the patient's charts. The liver, kidney, and muscle were analyzed with routine histology and electron microscopy. Mitochondrial function was assessed by respirometric analyses and enzymatic activity assays. Structure and sequence analyses of this specific variant were investigated by in silico analyses. Our patient showed the known features of the variant-specific phenotype, including macrosomia, congenital hyperinsulinism, transient hepatomegaly, and renal Fanconi syndrome. In addition to that, she showed liver cirrhosis, chronic kidney failure, and altered mitochondrial morphology and function. The clinical and biochemical phenotype had features of a new type of glycogen storage disease. Discussion: This case expands the p.(Arg85Trp) variant-specific phenotype. Possible pathomechanistic explanations for the documented multiorgan involvement and changes of symptoms and signs during development of this ultra-rare but instructive disorder are discussed.

Non-steroidal CYP17A1 Inhibitors: Discovery and Assessment.

CYP17A1 is an enzyme that plays a major role in steroidogenesis and is critically involved in the biosynthesis of steroid hormones. Therefore, it remains an attractive target in several serious hormone-dependent cancer diseases, such as prostate cancer and breast cancer. The medicinal chemistry community has been committed to the discovery and development of CYP17A1 inhibitors for many years, particularly for the treatment of castration-resistant prostate cancer. The current Perspective reflects upon the discovery and evaluation of non-steroidal CYP17A1 inhibitors from a medicinal chemistry angle. Emphasis is placed on the structural aspects of the target, key learnings from the presented chemotypes, and design guidelines for future inhibitors.

Loss of LGR4/GPR48 causes severe neonatal salt wasting due to disrupted WNT signaling altering adrenal zonation.

Disorders of isolated mineralocorticoid deficiency, which cause potentially life-threatening salt-wasting crisis early in life, have been associated with gene variants of aldosterone biosynthesis or resistance; however, in some patients no such variants are found. WNT/ß-catenin signaling is crucial for differentiation and maintenance of the aldosterone-producing adrenal zona glomerulosa (zG). Herein, we describe a highly consanguineous family with multiple perinatal deaths and infants presenting at birth with failure to thrive, severe salt-wasting crises associated with isolated hypoaldosteronism, nail anomalies, short stature, and deafness. Whole exome sequencing revealed a homozygous splice variant in the R-SPONDIN receptor LGR4 gene (c.618-1G>C) regulating WNT signaling. The resulting transcripts affected protein function and stability and resulted in loss of Wnt/ß-catenin signaling in vitro. The impact of LGR4 inactivation was analyzed by adrenal cortex-specific ablation of Lgr4, using Lgr4fl/fl mice mated with Sf1:Cre mice. Inactivation of Lgr4 within the adrenal cortex in the mouse model caused decreased WNT signaling, aberrant zonation with deficient zG, and reduced aldosterone production. Thus, human LGR4 mutations establish a direct link between LGR4 inactivation and decreased canonical WNT signaling, which results in abnormal zG differentiation and endocrine function. Therefore, variants in WNT signaling and its regulators should systematically be considered in familial hyperreninemic hypoaldosteronism.

Recognition of fold- and function-specific sites in the ligand-binding domain of the thyroid hormone receptor-like family.

Background: The thyroid hormone receptor-like (THR-like) family is the largest transcription factors family belonging to the nuclear receptor superfamily, which directly binds to DNA and regulates the gene expression and thereby controls various metabolic processes in a ligand-dependent manner. The THR-like family contains receptors THRs, RARs, VDR, PPARs, RORs, Rev-erbs, CAR, PXR, LXRs, and others. THR-like receptors are involved in many aspects of human health, including development, metabolism and homeostasis. Therefore, it is considered an important therapeutic target for various diseases such as osteoporosis, rickets, diabetes, etc. Methods: In this study, we have performed an extensive sequence and structure analysis of the ligand-binding domain (LBD) of the THR-like family spanning multiple taxa. We have use different computational tools (information-theoretic measures; relative entropy) to predict the key residues responsible for fold and functional specificity in the LBD of the THR-like family. The MSA of THR-like LBDs was further used as input in conservation studies and phylogenetic clustering studies. Results: Phylogenetic analysis of the LBD domain of THR-like proteins resulted in the clustering of eight subfamilies based on their sequence homology. The conservation analysis by relative entropy (RE) revealed that structurally important residues are conserved throughout the LBDs in the THR-like family. The multi-harmony conservation analysis further predicted specificity in determining residues in LBDs of THR-like subfamilies. Finally, fold and functional specificity determining residues (residues critical for ligand, DBD and coregulators binding) were mapped on the three-dimensional structure of thyroid hormone receptor protein. We then compiled a list of natural mutations in THR-like LBDs and mapped them along with fold and function-specific mutations. Some of the mutations were found to have a link with severe diseases like hypothyroidism, rickets, obesity, lipodystrophy, epilepsy, etc. Conclusion: Our study identifies fold and function-specific residues in THR-like LBDs. We believe that this study will be useful in exploring the role of these residues in the binding of different drugs, ligands, and protein-protein interaction among partner proteins. So this study might be helpful in the rational design of either ligands or receptors.

Variant predictions in congenital adrenal hyperplasia caused by mutations in CYP21A2.

CYP21A2 deficiency represents 95% of congenital adrenal hyperplasia (CAH) cases, a group of genetic disorders that affect steroid biosynthesis. The genetic and functional analysis provide critical tools to elucidate complex CAH cases. One of the most accessible tools to infer the pathogenicity of new variants is in silico prediction. Here, we analyzed the performance of in silico prediction tools to categorize missense single nucleotide variants (SNVs) of CYP21A2. SNVs of CYP21A2 characterized in vitro by functional assays were selected to assess the performance of online single and meta predictors. SNVs were tested separately or in combination with the related phenotype (severe or mild CAH form). In total, 103 SNVs of CYP21A2 (90 pathogenic and 13 neutral) were used to test the performance of 13 single-predictors and four meta-predictors. All SNVs associated with the severe phenotypes were well categorized by all tools, with an accuracy of between 0.69 (PredictSNP2) and 0.97 (CADD), and Matthews' correlation coefficient (MCC) between 0.49 (PoredicSNP2) and 0.90 (CADD). However, SNVs related to the mild phenotype had more variation, with the accuracy between 0.47 (S3Ds&GO and MAPP) and 0.88 (CADD), and MCC between 0.18 (MAPP) and 0.71 (CADD). From our analysis, we identified four predictors of CYP21A2 variant pathogenicity with good performance, CADD, ConSurf, DANN, and PolyPhen2. These results can be used for future analysis to infer the impact of uncharacterized SNVs in CYP21A2.

Loss of Protein Stability and Function Caused by P228L Variation in NADPH-Cytochrome P450 Reductase Linked to Lower Testosterone Levels.

Cytochrome P450 oxidoreductase (POR) is the redox partner of steroid and drug-metabolising cytochromes P450 located in the endoplasmic reticulum. Mutations in POR cause a broad range of metabolic disorders. The POR variant rs17853284 (P228L), identified by genome sequencing, has been linked to lower testosterone levels and reduced P450 activities. We expressed the POR wild type and the P228L variant in bacteria, purified the proteins, and performed protein stability and catalytic functional studies. Variant P228L affected the stability of the protein as evidenced by lower unfolding temperatures and higher sensitivity to urea denaturation. A significant decline in the rate of electron transfer to cytochrome c and thiazolyl blue tetrazolium (MTT) was observed with POR P228L, while activities of CYP3A4 were reduced by 25% and activities of CYP3A5 and CYP2C9 were reduced by more than 40% compared with WT POR. The 17,20 lyase activity of CYP17A1, responsible for the production of the main androgen precursor dehydroepiandrosterone, was reduced to 27% of WT in the presence of the P228L variant of POR. Based on in silico and in vitro studies, we predict that the change of proline to leucine may change the rigidity of the protein, causing conformational changes in POR, leading to altered electron transfer to redox partners. A single amino acid change can affect protein stability and cause a severe reduction in POR activity. Molecular characterisation of individual POR mutations is crucial for a better understanding of the impact on different redox partners of POR.

Synthesis and Structure-Activity Relationships of Novel Non-Steroidal CYP17A1 Inhibitors as Potential Prostate Cancer Agents.

Twenty new compounds, targeting CYP17A1, were synthesized, based on our previous work on a benzimidazole scaffold, and their biological activity evaluated. Inhibition of CYP17A1 is an important modality in the treatment of prostate cancer, which remains the most abundant cancer type in men. The biological assessment included CYP17A1 hydroxylase and lyase inhibition, CYP3A4 and P450 oxidoreductase (POR) inhibition, as well as antiproliferative activity in PC3 prostate cancer cells. The most potent compounds were selected for further analyses including in silico modeling. This combined effort resulted in a compound (comp 2, IC50 1.2 µM, in CYP17A1) with a potency comparable to abiraterone and selectivity towards the other targets tested. In addition, the data provided an understanding of the structure-activity relationship of this novel non-steroidal compound class.

Biased cytochrome P450-mediated metabolism via small-molecule ligands binding P450 oxidoreductase.

Metabolic control is mediated by the dynamic assemblies and function of multiple redox enzymes. A key element in these assemblies, the P450 oxidoreductase (POR), donates electrons and selectively activates numerous (>50 in humans and >300 in plants) cytochromes P450 (CYPs) controlling metabolism of drugs, steroids and xenobiotics in humans and natural product biosynthesis in plants. The mechanisms underlying POR-mediated CYP metabolism remain poorly understood and to date no ligand binding has been described to regulate the specificity of POR. Here, using a combination of computational modeling and functional assays, we identify ligands that dock on POR and bias its specificity towards CYP redox partners, across mammal and plant kingdom. Single molecule FRET studies reveal ligand binding to alter POR conformational sampling, which results in biased activation of metabolic cascades in whole cell assays. We propose the model of biased metabolism, a mechanism akin to biased signaling of GPCRs, where ligand binding on POR stabilizes different conformational states that are linked to distinct metabolic outcomes. Biased metabolism may allow designing pathway-specific therapeutics or personalized food suppressing undesired, disease-related, metabolic pathways.

Reducing FASN expression sensitizes acute myeloid leukemia cells to differentiation therapy.

Fatty acid synthase (FASN) is the only human lipogenic enzyme available for de novo fatty acid synthesis and is often highly expressed in cancer cells. We found that FASN mRNA levels were significantly higher in acute myeloid leukemia (AML) patients than in healthy granulocytes or CD34+ hematopoietic progenitors. Accordingly, FASN levels decreased during all-trans retinoic acid (ATRA)-mediated granulocytic differentiation of acute promyelocytic leukemia (APL) cells, partially via autophagic degradation. Furthermore, our data suggest that inhibition of FASN expression levels using RNAi or (-)-epigallocatechin-3-gallate (EGCG) accelerated the differentiation of APL cell lines and significantly re-sensitized ATRA refractory non-APL AML cells. FASN reduction promoted translocation of transcription factor EB (TFEB) to the nucleus, paralleled by activation of CLEAR network genes and lysosomal biogenesis. Together, our data demonstrate that inhibition of FASN expression in combination with ATRA treatment facilitates granulocytic differentiation of APL cells and may extend differentiation therapy to non-APL AML cells.

Characterization of Mutations Causing CYP21A2 Deficiency in Brazilian and Portuguese Populations.

Deficiency of 21-hydroxylase enzyme (CYP21A2) represents 90% of cases in congenital adrenal hyperplasia (CAH), an autosomal recessive disease caused by defects in cortisol biosynthesis. Computational prediction and functional studies are often the only way to classify variants to understand the links to disease-causing effects. Here we investigated the pathogenicity of uncharacterized variants in the CYP21A2 gene reported in Brazilian and Portuguese populations. Physicochemical alterations, residue conservation, and effect on protein structure were accessed by computational analysis. The enzymatic performance was obtained by functional assay with the wild-type and mutant CYP21A2 proteins expressed in HEK293 cells. Computational analysis showed that p.W202R, p.E352V, and p.R484L have severely impaired the protein structure, while p.P35L, p.L199P, and p.P433L have moderate effects. The p.W202R, p.E352V, p.P433L, and p.R484L variants showed residual 21OH activity consistent with the simple virilizing phenotype. The p.P35L and p.L199P variants showed partial 21OH efficiency associated with the non-classical phenotype. Additionally, p.W202R, p.E352V, and p.R484L also modified the protein expression level. We have determined how the selected CYP21A2 gene mutations affect the 21OH activity through structural and activity alteration contributing to the future diagnosis and management of CYP21A2 deficiency.

Finding New Molecular Targets of Familiar Natural Products Using In Silico Target Prediction.

Natural products comprise a rich reservoir for innovative drug leads and are a constant source of bioactive compounds. To find pharmacological targets for new or already known natural products using modern computer-aided methods is a current endeavor in drug discovery. Nature's treasures, however, could be used more effectively. Yet, reliable pipelines for the large-scale target prediction of natural products are still rare. We developed an in silico workflow consisting of four independent, stand-alone target prediction tools and evaluated its performance on dihydrochalcones (DHCs)-a well-known class of natural products. Thereby, we revealed four previously unreported protein targets for DHCs, namely 5-lipoxygenase, cyclooxygenase-1, 17ß-hydroxysteroid dehydrogenase 3, and aldo-keto reductase 1C3. Moreover, we provide a thorough strategy on how to perform computational target predictions and guidance on using the respective tools.

Characterization of Two Novel Variants of the Steroidogenic Acute Regulatory Protein Identified in a Girl with Classic Lipoid Congenital Adrenal Hyperplasia.

Congenital adrenal hyperplasia (CAH) consists of several autosomal recessive disorders that inhibit steroid biosynthesis. We describe a case report diagnosed with adrenal insufficiency due to low adrenal steroids and adrenocorticotropic hormone excess due to lack of cortisol negative feedback signaling to the pituary gland. Genetic work up revealed two missense variants, p.Thr204Arg and p.Leu260Arg in the STAR gene, inherited by both parents (non-consanguineous). The StAR protein supports CYP11A1 enzyme to cleave the side chain of cholesterol and synthesize pregnenolone which is metabolized to all steroid hormones. We used bioinformatics to predict the impact of the variants on StAR activity and then we performed functional tests to characterize the two novel variants. In a cell system we tested the ability of variants to support cholesterol conversion to pregnenolone and measured their mRNA and protein expression. For both variants, we observed loss of StAR function, reduced protein expression and categorized them as pathogenic variants according to guidelines of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology. These results fit the phenotype of the girl during diagnosis. This study characterizes two novel variants and expands the list of missense variants that cause CAH.

3ß-Isoobeticholic acid efficiently activates the farnesoid X receptor (FXR) due to its epimerization to 3α-epimer by hepatic metabolism.

Bile acids (BAs) are important signaling molecules acting via the farnesoid X nuclear receptor (FXR) and the membrane G protein-coupled bile acid receptor 1 (GPBAR1). Besides deconjugation of BAs, the oxidoreductive enzymes of colonic bacteria and hepatocytes enable the conversion of BAs into their epimers or dehydrogenated forms. Obeticholic acid (OCA) is the first-in-class BA-derived FXR agonist approved for the treatment of primary biliary cholangitis. Herein, a library of OCA derivatives, including 7-keto, 6-ethylidene derivatives and 3ß-epimers, was synthetized and investigated in terms of interactions with FXR and GPBAR1 in transaction assays and evaluated for FXR target genes expression in human hepatocytes and C57BL/6 mice. The derivatives were further subjected to cell-free analysis employing in silico molecular docking and a TR-FRET assay. The conversion of the 3ßhydroxy epimer and its pharmacokinetics in mice were studied using LC-MS. We found that only the 3ß-hydroxy epimer of OCA (3ß-isoOCA) possesses significant activity to FXR in hepatic cells and mice. However, in a cell-free assay, 3ß-isoOCA had about 9-times lower affinity to FXR than did OCA. We observed that 3ß-isoOCA readily epimerizes to OCA in hepatocytes and murine liver. This conversion was significantly inhibited by the hydroxy-Δ5-steroid dehydrogenase inhibitor trilostane. In addition, we found that 3,7-dehydroobeticholic acid is a potent GPBAR1 agonist. We conclude that 3ß-isoOCA significantly activates FXR due to its epimerization to the more active OCA by hepatic metabolism. Other modifications as well as epimerization on the C3/C7 positions and the introduction of 6-ethylidene in the CDCA scaffold abrogate FXR agonism and alleviate GPBAR1 activation.

Novel CYP19A1 Mutations Extend the Genotype-Phenotype Correlation and Reveal the Impact on Ovarian Function.

CONTEXT: The steroidogenic enzyme aromatase (CYP19A1) is required for estrogen biosynthesis from androgen precursors in the ovary and extragonadal tissues. The role of aromatase, and thus estrogens, is best illustrated by genetic variations of the CYP19A1 gene leading to aromatase deficiency or excess. OBJECTIVE: The objective of this work is to characterize novel CYP19A1 variants. DESIGN SETTING AND PATIENTS: Variants causing aromatase deficiency were suspected in four 46,XX children of African and Indian origin by careful clinical phenotyping. Sequencing of the CYP19A1 gene identified novel variants. Minigene experiments, aromatase activity assay, and computational, and histological analysis were used to characterize the variants. MAIN OUTCOME MEASURE AND RESULTS: CYP19A1 variants were found in all patients: a deletion in intron 9 leading to p.P423_H503del, a delins variant at p.P154, and point variants p.V161D, p.R264C, p.R375C. Except for R264C, all variants showed a loss of function. Protein structure and dynamics studies were in line with functional assays. The 2 female patients with delins variants manifested with ambiguous genitalia at birth. Histologic investigation revealed normal ovarian tissue on one side and a streak gonad on the other. Two female patients presented with abnormal pubertal development and polycystic ovaries. CONCLUSION: In girls, aromatase deficiency usually manifests at birth, but diagnosis may also be made because of abnormal pubertal development or ovarian torsion due to (poly)cystic ovaries. The ovary harboring CYP19A1 variants may present as streak gonad or appears normal at birth, but is then at very high risk to produce cysts with aging and is therefore prone to ovarian torsion.